R=660;
L=100E-6;
tau=L/R;
dt=tau/100;
t=0:dt:6*tau;
u1(1)=0;
u1(2)=0;
u2(1)=0;
lt=length(t)
u1=randn(1,lt);
for count=2:1:length(t)

u2(count)=(u1(count)+(tau/dt)*u2(count-1))/(1+tau/dt);    
end;
plot(t,u2/3)